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Abstract. The energy minimization problem associated to uniform, isotropic, linearly elastic rods leads to 
a geometric variational problem for the rod centerline, whose solutions include closed, knotted curves. We 
give a complete description of the space of closed and quasiperiodic solutions. The quasiperiodic curves are 
parametrized by a two-dimensional disc. The closed curves arise as a countable collection of one-parameter 
families, connecting the m-fold covered circle to the n-fold covered circle for any m, n relatively prime. 
Each family contains exactly one self-intersecting curve, one elastic curve, and one closed curve of constant 
torsion. Two torus knot types are represented in each family, and all torus knots are represented by elastic 
rod centerlines. 



Introduction 

The elastic curve of Bernoulli and Euler arises as the solution of the variational problem of minimizing 
the total squared (geodesic) curvature (defined below) of a curve with prescribed boundary conditions and 
fixed length. In [pj, Langer and Singer parametrized the space of (similarity classes of) elastic curves by 
a triangular region of the plane; the closed curves form a countable set of points lying on a single curve c 
within the triangle. The elastic curves determined by points on c are quasiperiodic; that is, they are either 
closed curves or they wind densely around a torus of revolution and are self-congruent under a discrete group 
of rotations of the torus about its axis. The closed curves are in 1-to-l correspondence with the knot types 
which are {k, n)-torus knots for k < n/2. 

A generalization of the elastic curve can be obtained by adding the additional constraint that the integral 
of the torsion be held constant. Thus we seek an extremal 7 for the Lagrangian 

T[-f] = Xi ds + X2 Tds + Xs K^ds, A3 ^ (1) 
J y J ^ J ^ 

Here, k and t, the curvature and torsion, are the coefficients of the Frenet equations for a canonical or- 
thonormal (Frenet) frame T, N, B along 7, given by 

^ ^T, T' ^ kN, N' = -kT + tB, B' = -tN. 
ds 

The constants Ai, A2 and A3 are Lagrange multipliers. 

That this is an appropriate generalization of elastic curves is a consequence of a theorem of Langer and 
Singer (cf. Theorem 1): If 7 is the centerline of a uniform symmetric Kirchhoff elastic rod (see section 
^), then 7 is an extremum for the Lagrangian In fact, as we show in section the converse of this 
theorem also holds. Thus we may refer to extremals of (0) as elastic rod centerlines. Our main results are 

Theorem 1. Every torus knot type is realized by a smooth closed elastic rod centerline. 



Theorem 2. The similarity classes of quasiperiodic elastic rod centerlines are parametrized by the closed 
unit disc. (Quasiperiodic centerlines are those which are either closed or wind densely around a torus.) The 
closed rod centerlines form a countable family of curves in the disc. Each such curve represents a regular 
homotopy of closed curves within the set of smooth closed elastic rod centerlines. For any relatively prime 
positive integers k, n, there exists a regular homotopy, between the fc-times-covered circle and the n-times 
covered circle. The homotopy includes exactly one elastic curve, one self-intersecting elastic rod centerline, 
and one closed curve of constant torsion. 

Figure 1 shows the members of one such homotopy of closed elastic rod centerlines. 
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Figure 1. Elastic rod centerlines along the liomotopy joining the once-covered circle to the 
twice-covered circle. 



We will now give an outline of the rest of the paper. In §1 we detail the explicit solution of the variational 
problem for (|l|). In §2 we show how parameters can be chosen so that the equilibria are smooth closed 
curves, and we show that, with a few exceptions, these are embedded and lie on tori of revolution. In §3 
we introduce a new parametrization of the space T) of quasiperiodic elastic rods, and point out loci inside 
this space that represent elastic curves, curves of constant torsion, self-intersecting curves, etc. In §4 we 
consider what the elastic rods look like at the boundary of T), and prove Theorem ^ In §5 we examine the 
level curves of a certain function on T> which turn out to give the homotopies of Theorem ^, and prove a 
uniqueness theorem (Theorem^, conjectured in [Q, on the knot types of elastic rods of constant torsion. In 
§6 we relate the variational problem to the Kirchhoff elastic rod and discuss the question of stability. The 
problem of stability for closed elastic rods has been extensively investigated, using bifurcation theory, by K. 
Rogers (see H]). Here we examine the second variation formula in certain specific examples. This leads 
to a result in §7, specifying the precise nature of stability for the untwisted "figure eight" elastic rod. 

The authors would like to thank Joel Langer, Annalisa Calini, and Ron Perline for helpful discussions 
and encouragement. 

1. Elastic Rods in Cylindrical Coordinates 

In this section, we will recapitulate the integration of the Euler-Lagrange equations for the Kirchhoff 
elastic rod by Langer and Singer ||^. 

Suppose 7 is a space curve which is critical for the functional (|l|) with respect to variations that keep the 
endpoints fixed to second order. The resulting Euler-Lagrange equations (see §4 in |P|) have the following 
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first integrals: 

K2(2A3r- A2) =c 

The latter equation allows us to express the curvature as 

2 

2 2ri P 2u W 4- '^0'^ 

K = Ko 1 2™ ^'W)' 

2w 

where < p < w < 1, and p, w are related to the other parameters by 

4A1A3 - A2 /Q 2 2 IN 

= — W(3w - p - 1) 



(2) 



A§ 



(3) 



|Jp^M^ = ±(A3.^-2Aif + ^(^-A2f (4) 



Another consequence of the first integrals is that the vector field 

J = i(A3K2 - 2Ai)r + X^K.N + K(A3r - A2)S 
is constant along 7 (cf. Thni. 2 in ||^). Setting k = in (||) gives 

I2 

In addition to J, 

/ = X2T + X3KB 

is the restriction to 7 of a Killing field, generating a screw motion in R'^. Since / + J x 7 is constant, J turns 
out to be parallel to the axis of the screw motion. In fact, we will translate the curve so that I = aJ + j x J. 

The Killing fields / and J define a natural system of cylindrical coordinates r, 0, z, in which J points along 
the positive z-axis. Computing / • J in two ways gives 

a/i^ = iAsc - A1A2. 

Computing / • / in two ways gives 

^l'{r' + a')^Xl + \lK', (5) 
giving us the r-coordinate along 7. Computing T ■ J gives 

= ^X3k'^ - Ai (6) 

and computing T ■ I gives 

a^izs - A2 
Os = 2 • 

Since and are both of the form A + B sn^ t, z{s) and 9{s) can be computed in terms of elliptic integrals 
of the second and third kind, respectively; we will do so below. 

2. Closure Conditions 

Without loss of generality, we may assume A3 = 1 . By applying a similarity to a rod we may assume that 
the maximum curvature Kq = 1. Then the shape of the elastic rods depends on Ai, A2 and the constant of 
integration c, with /i > determined by (^). (We will ultimately use a different set of parameters, though.) 
We would like to see what choices of the parameters lead to closed curves. 

Since is a 2if-periodic function of t = kos/{2w) (where K = K{p) is a complete elliptic integral of 
the first kind), equation ^) implies that r is periodic. Likewise, (^) shows that zt is periodic, and z will be 
periodic if and only if ztdt ~ 0; this is equivalent to 

2 2 I 2\2 2i?(p) 

(cf. equation (26) in [^). So, in order for a rod to have periodic z coordinate, 
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must not be negative. Since the behaviour of A{p) as a function of p wih be important in what foUows — for 
example, the sign restriction gives us an upper bound on p — we will summarize the properties of A{p) which 
we will need. 

Proposition 2.1. For p £ [0, 1], A{p) is a smooth monotone decreasing function of p, with A{0) — 1 and 
A(l) = -1; and, 

iovpe (0,1). 

Proof. The behaviour of A at each end follows from that of E and K. Using differentiation formulas from 
m, we have 



dA 
dp 



E-K 
pK 



and 



Finally, 



E 



dp 



E-{p'YK 

p{p'Y 



{A + p^) 



2 

P 

E 
'Kj/ 



{1~E/Kf + 



< 0. 



p^E^ 



{p')^K^ 



< 



2,-(A + 2/-l) 



V-(A + 2p^-l)7(p') 



l\2 



The last equation shows that, after A + 2p^ — 1 = at p = 0, it remains positive until p' — 0. 



□ 



Let Pniax ~ .9089085 be the value of p where ^ = 0. When p — Pmax, the closure condition (|^) forces 
w = p, hence c = and the torsion is zero, and A2 = 0; then the elastic rod becomes a planar figure-eight 
elastic curve (cf. Figure 1 in [||.) 

When we impose the closure condition (|^) , the formula for the derivative of 9 becomes 

(«) 



dt 



1 - M sn2 i ' 



where 



U = a^w, 



N ■ 



^{l~A)-2\2tiw^ 



p^ + A ~ p^ 

(It will follow from calculations in the next section that < A — < {p'Y , so p^ < M < 1.) Assuming 
M < 1, we now use formula 434.01 in to compute 9. The following identity is necessary for simplifying 
the result: 



M 



P 



A-U^' 



(1 -Af)(M-p2) 
M ■ 



(9) 



Then 



where 



9{t) ^Ut + sgniN) 



E.^l]F{^^p') + E{^,p')]t^ir\n 



K 



2i 



e,{t-tF{lp')) 
ei{t + iFilp')) 



^ = sin ^ 



{p'^M 



(10) 



p' — and 81 is the Jacobi theta function with period 2K and Qi{K + iK') = 0. The change in 

9 over one period of the other coordinates r, z is given by 

(11) 



A9 = 2KU + 2 sgn(iV) 1^{E - K)F{^,p') + KEi^, p') j 
The rod will be periodic if the parameters are chosen so that A9 is a rational multiple of 27r. 

Theorem 3. An elastic rod center line of non-constant curvature can intersect itself only at the origin of 
the natural system of cylindrical coordinates. If this does not occur, and the centerline is closed, then it is 
embedded and lies on an embedded torus of revolution. 
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Proof. Note that the curvature is constant if and only if p = 0, giving a planar curve; we will assume p > 0. 
We've observed that and z are 2if -periodic functions of t. The minimum value of is non-negative (see 
Lemma 3.1 below) and occurs when t = 0. 

Suppose that rmin > 0. Then r{t) is a smooth positive function with stationary points at t = and t — K, 
giving the maximum and minimum values in each 2K period. Because of the condition (|7|), 

where 9 is the Jacobi theta function with period 2K and Q{iK') = 0. Then z{t) is a smooth odd function 
with stationary points where dn t = E/K, giving the maximum and minimum values in each period. Thus, 
r{t) and z{t) trace out a simple closed curve in the rz half-plane. It follows that the rods lie on embedded 
tori of revolution. 

We would like to show that the rod is in fact embedded. Suppose 6 and (f) are longitudinal and latitude on 
the embedded torus. We have already seen that 4>{t) is a monotone increasing function along the rod, with 
(j){t + 2K) = (t){t) -\- 2tt. Since Ot is periodic, it follows that the rod can intersect itself only when it closes up 
smoothly. 

In the case where rmin = 0, we can similarly argue that r{t) and z{t) trace out a closed curve in the rz 
half-plane, which is smooth except at the origin. Replacing cf) with a local coordinate along the smooth part 
of this curve, we again see that the rod cannot intersect itself except at the origin. □ 

We conclude this section by giving formulae for the rods in Cartesian coordinates. Let A = A6/{2K). 
Then 



l2Kpp' e'^^QAt ± iF) + e-'^*ei{t T iF) 

r cos — * ' 



2p,wQ{t)H,{iF) ^^^^ 



l2Kpp' e'^'Qiit ± iF) - e-'^'SiitT iF) 

rsmy = -\/ — , 

71- 2i^iwQit)HiiiF) 

where F = F{^,p'), and the upper sign is taken when A'' > 0, the lower sign when N < 0. 

The formula ( p^ ) for z bears a striking resemblance to the representation of planar elastic curves in terms 
of theta functions obtained by Mumford In fact, one can check that as p Pmax, and our elastic rod 
becomes planar, ( p^ ) and ( p^ agree with Mumford's formula in the limit. 

3. The Space V of Quasiperiodic Elastic Rods 

Because the curvature and torsion are 2if -periodic in t, the curve will consist of successive congruent 
segments which join together smoothly; when we impose the closure condition in the z coordinate, the 
segments will be congruent by a rotation in 9. Thus, the rod centerline will either close up smoothly after a 
finite number of segments, or wind densely around a torus. We call these quasiperiodic elastic rods. 

Such curves depend on two parameters; one possible choice of parameters is the elliptic modulus p and 
the parameter A2. However, we will view the parameter space as a unit disk in the XY plane, where X and 
Y are related to the other parameters as follows: set 

2 

x = X2w, y= =, z = cw^. 

Then, from (|) and (g), 

w = y/Y^+p^, Z = Yy/{p')^ ~Y^, X^ + Y^ = A{p). 

(The last equation is equivalent to ^.) If p £ (0,pniax) then w > 0, and X and Y are seen to be smooth 
invertible functions of c and p, tending to the origin as p approaches Pmax and tending to the edge of the 
unit disk as p approaches zero. In terms of X, Y, Z and w, the constants used in the rod formulae are: 



1 A-l A)^ +'i(Z -wX)"^ 



2 4w2 ' 4u)2 

2w{Z-wX)+X{l-A) 



U = apw 
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We will let T) stand for the open unit disk in the XY plane, with the origin omitted for convenience. On 
2?, we will also use the angle coordinate 0, where 

X = \/A{p) cos Y = ^/A{p) sin (j). 

Thus, p, (j) are also valid coordinates on T). 

There are several interesting loci inside T>. First, X — gives elastic curves. Next, Y — gives elastic 
rods of constant torsion r = A2/2; these rods are the starting point for the more elaborate knots of constant 
torsion produced in Next, recall that the Killing field 

/ = A2T + X3KB = aJ + J X J 

generates the motion of the elastic rod, up to a tangential piece, under the Localized Induction Equation (cf. 
1^, Theorem 3 and Corollary 5.1). While in general this is a screw motion, it is natural to ask when it is a 
pure rotation. Since J is the translational Killing field, defining the z axis for the cylindrical coordinates, / 
generates a rotation when / • J = 0, i.e. when a = 0. Setting a = gives 

2Zw = X{A - 1 + 2p2 + 2Y^), 

showing that X and Y must have the same sign along this curve. Substituting in polar coordinates on the 
disk gives 

. 2 (A-l + 2pY 

'P i_A2-2A(A-l + 2p2)- 

Kida [Q was the first to investigate space curves that move by a pure rotation under the LIE, and for 
this reason we call this curve in V the "Kida curve". It is amusing exercise to verify that, as p ^ 0, sin^ 
approaches 2/3 along the Kida curve (see Fig. 2). Kida also claimed, without proof, that the closed elastic 
rods in this family have the knot type of a (fc, n) torus knot with \n/k\ > 2. Our results below confirm this. 



Figure 2. The parameter space V of quasi-periodic elastic rods. Points along the Y axis 
correspond to elastic curves, along the X axis to elastic rods of constant torsion, and along 
the S-shaped curve to self-intersecting rods. The Kida curve runs through the first and third 
quadrants. 

A fourth curve in V of geometric interest corresponds to those rods which are self-intersecting. To locate 
this curve in I?, we need the following 

Lemma 3.1. The minimum value of along an elastic rod is given by 

rL, = (A-C/2)/(^W) 

If ^ is the angle in (^^, then 

yl-J72 = (l-p2)cos2| = y2 (14) 
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where 

V := {2Xw^/l-w^ -Y{l + A- 2v?)) / {'i^iv?) . 

Proof. Since 

r'^^l^u? ^X^ + II?- -U"^ ST? t, 

the minimum value of the right-hand side is A — [/^. The proof of ( |l4| ) is a simple (but tedious) calculation 
using the equations between various parameters, and we will omit it here. □ 

From Q, we see that r = is possible only when V ~ 0, i.e. 



2Xw\/l~w^ = y(l + A - 2p^ - 2Y^) 
Substituting in the polar coordinates on V and squaring gives 

cot0 = ± . 15 

2pp' 

Substituting back into the equation for X and Y shows that we should choose the minus sign in ( p^ , giving 
a curve in the second and fourth quadrants of T). One can check that, in Figure |[ F is positive to the right 
of this curve and negative to the left. 

It will be shown in §4 that V"^ < [p'Y on T). If we set cos^ ~ V/p' , then angle ^ is (up to 2ti) an analytic 
function of X and Y that is equal to ^ to right of the V — curve and equal to tt — <^ to the left. (Remember, 
^ only has range [0,7r/2].) It follows from and (|lj) that 



y2 + p2 ■ 

Thus, = if and only \iV — 0. In fact, the sign of N is opposite that of V . (One can check this, say, by 
evaluating N along the Kida curve, where [/ = 0.) Since the term 

2({E-K)F{^,p')+KE{i,p'] 



in (11) has limit tt as ^ — > 7r/2, we see that 

M = 2KU - 2{E - K)F{i,p') - 2KE{^,p') (16) 
plus an additional 27r when cos^ < (see Figure 3). 




Figure 3. A plot of A6'/(27r) as a function of (p, for p — 0.8. Note the discontinuities, 
corresponding to rod centerlines that are self-intersecting. 
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4. Behaviour on the boundary of V 

As we approach the boundary of I?, the eUiptic modulus approaches zero, and the equation for the 
curvature 

= 1 - sn^ t 

shows that the curvature of the rod becomes constant, provided w has a nonzero Umit at the boundary point. 
(Henceforth we will assume this is the case; at the points where p = w = 0, curvature has an undefined 
limit.) The torsion has limiting value 

r = (A2 + c/«2)/(2A3) = (Z + wX)lw^ 

But along the boundary, 

Z -^Y\/l-Y^, wX ^ X\Yl 

so the limiting value of Z + wX is Y\X \ + X\Y\. Thus, in the second and fourth quadrants the torsion goes 
to zero, and the rod is a circular arc. In the first and third quadrants the torsion has a nonzero limit, and 
one might conclude that the limiting shape of the rod that of a helix. However, the rod itself goes off to 
infinity in the cylindrical coordinates, i.e. 

One can see this by checking that 

lim(l - = 1, lim2(Z- wX) = -cot(/) 

and hence [/ — > and /i ^ along this part of the boundary. What one observes experimentally is that, as 
the parameters values approach the boundary in the first and third quadrants, the congruent segments of 
the rod become shorter (as A6' approaches zero) and more and more helical (see Figure 



Figure 4. Three congruent segments of a nearly helical elastic rod centerline (p = 0.001, 
cj) = 7r/4). 

For what follows, it will be necessary to know how behaves along the boundary of V. It is easy to 
verify that, as we approach the part of the boundary that is strictly inside the second or fourth quadrants, 

C/^y, ^lw'^-^\XYl M^O, N ^Y 

and so 9t 2Y. Hence along this part of the boundary, 

AO = 2nY. (17) 

In particular, as (f) increases, A9 decreases from 2tt to along the second quadrant edge, and increases from 
— 27r to along the fourth quadrant edge. 

As we approach the part of the boundary strictly inside the first or third quadrants, 

Ai^O, C/^0, N^O 

and hence 6t ^ 0. 
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Proof, [of Theorem |l[ Because of ( p^ ) and the continuity of A9, we see that for any rational number m/n 
such that \m/n\ € (0, 1), there are points in V, near the edge in the second (m/n > 0) or fourth (m/n < 0) 
quadrants, where AO = 2iTm/n. By Theorem |[ the corresponding elastic rod will be embedded and lie on 
a torus of revolution, and will complete m circuits around the 2;-axis and n circuits around the waist of the 
torus before it closes up smoothly. □ 

5. Behaviour of AO on the interior of V 



Let AO stand for the smooth part of AO (i.e. the right-hand side of (|16[)). In this section we will show that 
AO has smooth level curves in T) that cross the disk from the second quadrant edge to the fourth quadrant 
edge, and which cross each of the X-axis, the F-axis, the Kida curve (where J7 = 0), and the self-intersection 
locus (where ^ = 0) exactly once. We will start by showing that p and U can be used as coordinates over 
most of T), and establishing the smoothness of the level curves in those coordinates. 

From (p^), we have A~ — {p'Y cos^ ^, and hence for a fixed p, 

UdU/d(t>^{p' f cos ^ sin ^d£,/d(t>. (18) 
Lemma 5.1. For p e (0,pinax) fixed, dU/dcj) = if and only if cos^ = 0. 

Proof. Note that V is the numerator of cos^, and, on the interior V, V =^ ii and only if cos^ = 0. Using 
the derivative formulae 

dX/d4>^~Y, dY/d(t)^X, dw/dcj) = XY/w, 

we obtain 



dU _ V{2XZ - 2wX^ + w{A - l)){2Y{wX - Z) + ^l-ui^{A - 1)) 

Note that V is the numerator of cos^, and, on the interior T), V = ii and only if cos^ — 0. 

To complete the proof, we need to check that the other two factors in the numerator of dU /d(j) cannot be 
zero. The middle factor is zero if and only if 



squaring both sides gives 



2XY\/l~w^ = w(l + X^- Y^); 
w'^{l-Ay +4:p^X^ = 0, 



which is impossible for p G (0,pmax)- The last factor in (19) is zero if and only if 

2wXY = Vl - u;2(l - A + 2Y\ 

squaring both sides gives 

(1 - 'u?){\ - Af + 4^2(1 - ^ - /) 0, 
which is also impossible for p £ (0,pmax), because of Prop. □ 

Corollary 5.2. For p e (0,pmax), dcos^/dcj) = if and only if [/ = 0. 

Proof. Since cos^ = V/p', putting ( |l9|) into (|l8|) allows us to cancel out V on both sides. Then the result 
follows from the remaining equation. □ 

Corollary 5.3. For p e (0,pinax), sin^ ^ 0. 



Proof. If sin<^ = 0, then by (|18|), either U — or dU/dcj) = 0. The previous lemma rules out the latter, and 
we just have to verify that, along the Kida curve, cos^ ^ 7^ 1- 
When U = 0, 2wZ = X{A - 1 + 2w^), and hence 

^ A{A~1 + 2p2) 
Ap'^w^Y 

and 

16Ai2w4 = (1 _ Af + 4A(1 - A - _p2). 
But squaring each side of 2wZ = X{A — 1 + 2w'^) gives 

A{A-l + 2p2)2 = Y^iiA - 1 + 2pY + 4(A + p'^){l -A- p^)) = IGfj.'^Y^w'^. 
Hence cos^ ^ = A/{p' f < 1. □ 
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By Lemma 5.1, we know that away from the curve V = 0, p and U can be used as coordinates on V. 
Using (Jig) and (p')^ cos^ — A — U"^, we compute 



1 dAe 
2K~dU 



1 



U{ElK-{p'Ysin' 



[p'Y cos ^ sin ^i / 1 — [p'Y s\T? ^ 



(20) 



Lemma 5.4. If p G (0,pmax) is fixed and cos^ ^ 0, d/S.9/dU = only when X = 0. 

Proof. Let U = ^/Acoslo and V — p'cos^ = V^sinoj. Using (p')^sin^^ = {p'Y — ^sin^w, the right-hand 
side of (E^) is zero if and only if 



tsaiujJ{{p'Y - Asin^w)(p2 + Asin^cj) - {E / K - {p'Y + Asin^ uj) = 0. 



(21) 

We will show that the left-hand side of (|21| ) is a monotone function of uj. Then it will follow that the 
right-hand side of (po|), as a function of angle 0, can be equal to zero at most once between its vertical 
asymptotes, which occur when cos^ = 0. 

That (EQ) vanishes when X = is easy to verify. Let 



Then, when X = 0, = A, and 



tan w 



-{l + A-2w^) 
2w\/T 



h — p^ + Asin^ u). 



{l-Af 



I6p?w^ ' 



1 - h 



{l-w^){l + Af 
l&p?w^ 



and 



IQp^w^ = (1 - Af + 4y2(i _ 

Substituting these in the left hand side of (|l]) gives zero. 

To complete the proof, let / ~ i&nuj^hil — h) — h, which differs from ( |2l| ) by a term depending only on 
p. We compute 



dj h{l ~ h) sec^ Lo + A{1 - 2h) sin^ w - 2Asinwcosu;^/i(l - h) 
'dui ^ ^h{l - h) 



_ {\Jh{l - h) - AsmujcosLvf + tan^ lj{A + p'^){l - A-p^) 

Since all the terms in the numerator and denominator arc nonnegative, and never simultaneously zero, we 
are done. □ 



Lemma 5.5. For p E (0,pmax) fixed, A9 is a monotone increasing function of (f> when X > 0, and monotone 
decreasing when X < 0. 

Proof. First, consider just the X > side. Suppose we are at a point where cos^ ^ 0; then by ( pTl| ) and 
( ^.4| ), dA9/d(p 7^ 0. On the other hand, when cos^ = 0, this happens in the fourth quadrant, and by (^^), 
d^/d<j) ^ there. From jl^ ) we have 



1 dAQ dU di 



2K 



(E/K-l) 



'1- (pO^sin^e 
f A + 2p^ - r 



1 - [p' Y SIT? £, 



2p 



^0. 



To get the sign for dAO/dcj), we check that when [7 = 0, AO is an increasing function of J7, and, by checking 
the signs of each term in (|l9|), that U is an increasing function of (j) when cos^ > 0. 

For the X < Q side, the argument is similar, except that we find that ?7 is a decreasing function of when 
[7 = 0. □ 
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Proof, [of Theorem ^ Given k and n coprime, we will show that there exists a curve in V along which 
A6 — — 27rfc/(fc + n), running from the fourth quadrant of P to the second quadrant (see Figure ^). This 
curve will intersect the V = curve once, afterwhich AO = 27rn/(fc + n). At each point along the curve, we 
form a closed elastic rod by letting t run through k + n periods of length 2K. (This is the least number at 
which A6 is an integer multiple of 2tt.) Thus at one end the rods will converge to an fc-covered circle, and 
at the other end to an n-covered circle. 

Assume first that k/{k + n) < 1/2. For elastic curves (corresponding to the line X ^ Q in T)), it was 
shown by Langer and Singer that AO is a monotone function of p, decreasing from to — tt as we approach 
the origin along the positive F-axis (cf. Figure 2 and Appendix). We will follow the level curve along 



which A9 = — 27rfc/(fc + n) from a point along the positive Y-axis. Because of (5J_) and (5^), the level curve 



is perpendicular to the Y-axis; however, because of (5.5), as we follow it into the first quadrant, p must 



be strictly decreasing along the curve. It cannot veer toward the edge of the first quadrant, since AO has 
limiting value zero there. Thus, it crosses down into the fourth quadrant, as p continues to decrease. 

In the fourth quadrant, the level curve cannot cross the curve V = 0. To see why, note that ^ = gives 
^ = 7r/2 and so 

Ae\v=o= '^KU -2{{E- K)K' + KE') = 2KU - vr, 

using Legendre's relation. Since C/ < in the fourth quadrant, AO < —tt along the V = curve. Hence the 
level curve must end at the appropriate point on the edge of the fourth quadrant. 

As we follow this level curve of AO in the other direction, from the positive F-axis into the second 
quadrant, we see that p must again be monotone decreasing. If it continues into the third quadrant, it 
cannot meet the third quadrant edge because AO = there; if it continues into the fourth quadrant, p must 
now be increasing, and the curve is forced into the origin since it cannot cross the V = curve. This too is 
impossible, since AO = — tt at the origin. Thus, the curve must end at the appropriate point on the edge of 
the second quadrant. 

It is clear that this level curve must cross the curve V — at least once in the second quadrant. Since 
(|l|) shows that = A when 0, 



AO\v=o= ±2^/K{2E - K) - tt. 

Then, since K(2E — K) is a monotone decreasing function of p, AO is a monotone function along the curve 
V = Q, and the level curve crosses V = Q only once. 

If k/{k + n) > 1/2, then we can find a level curve along which AO = —2Tm/{k + n), passing from the fourth 
quadrant, through the first, to a point in the second quadrant where AO — 2tt — 2?™/ (k + n) — 27rfc/(fc + n). 
Rotating this level curve by 180 degrees about the origin, we obtain a level curve of AO, passing from a point 
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in the fourth quadrant where = — 27rfc/(fc + n), through the third quadrant, to a point in the second 
quadrant where = 27rn/(fc + n). (Recall that Ad is an odd function on V.) Once again, this curve 
intersects the curve V = exactly once. 

That these homotopies contain exactly one elastic curve and one curve of constant torsion follows from 
the monotonicity of A0 along the X = line, and Theorem ^ below. □ 

The families of elastic rods in obtained in the proof of Theorem || all have the same discrete symmetry 
group, generated by rotation by 27rfc/(fc+n). Among these rods, there is one elastic curve, one rod of constant 
torsion, one self-intersecting centerline, and one curve of Kida type moving by a rotation under the LIE. The 
elliptic modulus is greatest for the elastic curve; since the length, with respect to the 'natural' parameter t, 
is 2{n + m)K, one might say that for these rods of fixed symmetry and fixed maximum curvature kq = 1, 
the 'natural length' is greatest for the elastic curve. This seems to echo the mountain-pass argument in 
0, which showed that among curves with fixed symmetry and fixed length, this same elastic curve is the 
minimax for J K^ds. 

Theorem 4. Given any relatively prime integers k,n such that \k/n\ < 1/2, there exists a unique smooth 
closed elastic rod of constant torsion with the knot type of a (k, n) torus knot. 

Proof. Existence follows from the continuity of A9 along the X-axis in P, and its limiting values of — tt as 
X ^ 0^ and n as X ^ 0^ . To get uniqueness, we need to verify that A9 is a monotone function of p along 
the X-axis. Lemma 5.1 and ( [l^ ) imply that p, ^ are valid coordinates in the vicinity of this axis. We will 



differentiate A9 with respect to p, holding ^ fixed. 

From Q, using the differentiation formulas from the proof of Prop. |2.1| 

dU 2 / \ 

Since can be written as A6 — 2KU — 7rAo(^,j3), where Aq is Heuman's lambda function (see Q, 150.03), 
then 

IdAe _ ( E K\ K ( E^ 2 2 A (K - E)sm^cos^ 

(K-E) I (p')2 cos2 ^ - 1 - C/2 sin ^ cos ^ 



If this is zero, then 



([/2 + 1 - {p'f cos2 e)(l - (p')^ sin2 e) = 2C/(p')2 sin^cosCy^l- (pO^sin^^, 
which is equivalent to the vanishing of a sum of squares: 



2 



Uyi - (pO^sin-'C- (p ) sin^cosCj +p=0. 

Of course, this does not vanish on V. We have thus shown that, away from the U = and V = curves, 
dA9 /dp 7^ 0. We remark that this monotonicity calculation includes that of as a special case. □ 



6. Elastic Energy and Stability 

In the rest of the paper, we will discuss in more detail the relationship between critical curves for the 
geometric functional T and the centerlines of elastic rods. Of particular interest to us will be the question 
of which, if any, of the closed nonplanar critical curves for T are centerlines for stable elastic rods, i.e., local 
minima for the appropriate energy functional and boundary conditions. 

One can model the stress experienced by a rod composed of some elastic material by associating to the 
configuration of the rod a material frame (Mi, M2, T). This is an oriented orthonormal frame such that T is 
tangent to the centerline of the rod and Mi and M2 track the twisting of the material about the centerline. 
The components Ui of the Darboux vector for this frame, defined by 



T' ^ -U2M1 + M1M2, M[ = U2T - -U3M2, = -uiT + U3M1, 
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are called the material strains. (The prime denotes derivative with respect to arclength.) In the standard 
linear model, the stress experienced by a uniform, isotropic elastic rod is given by the energy 

£= a{ul + ul) + fSujds ^ / aK^ + Pujds, (22) 
Jo Jo 

where L is the length of the rod and a and /3 are nonnegative constants depending on the physical character- 
istics of the rod. Equilibria for £ are known as Kirchhoff elastic rods. When (3 — as well, £ only depends 
on the curvature of the centerline of the rod, and equilibrium centerlines are elastic curves. Equilibrium 
equations for ( p^ ) are classical (see |8| for a modern derivation) , and imply that the twist rate is constant 
along equilibria. 

Now we will formulate a boundary value problem for £ that is appropriate to rods with closed centerlines. 
Suppose the centerline 7 is a curve which possesses a smooth generalized Frenet frame (i.e., satisfying the 
usual Frenet equations and extending smoothly across inflection points). Let ip be the angle such that 
Ml = N cosip + Bsmip. Then the total change Aip of this angle, measured by integrating -0', measures the 
twisting of the material frame relative to the geometry of the centerline. It is related to the twist rate M3 by 
the balance equation: 

U3ds^Atp + Jrds. (23) 

The boundary conditions we will consider are requiring the centerline of the rod to be closed, and prescribing 
the value of Aip. (Physically, this can be achieved by attaching collars to the ends of the rod, joining the 
ends together, and then adjusting the total twist by rotating one collar against the other.) The following 
theorem clarifies the relationship between critical rods for this boundary value problem and critical curves 
for §). 

Theorem 5. Given any constants a,f3 > 0, and any curve 7 which is critical for a geometric functional 
of the form (|l|), there is a rod with centerline 7 and twist rate M3 equal to a constant m, which is an 
equilibrium for £ with respect to variations that preserve a certain prescribed value for Aip. The condition 
a : 2/3m = A3 : A2 determines m, and the value for Aip is determined by (|2^). 

Proof. We will at first assume that 7 has no inflection points. Let L be the length of 7, and let TZ be the 
set of rods whose centerlines are curves of length L with no inflection points, and where Aip has some flxed 
value. On 7, construct a rod F with this prescribed value and M3 — m, determined by 



Lm = A-0 + J Tds (24) 
Let 5 stand for the gradient operator for functionals on TZ. At F G 72., 

6 J (1/3)^ ds — 2m 5 j u^ds — 2m ^ j '^s, 
the last equality obtained by differentiating (p3|); hence, 

S £ — a S ds + 2(3m 5 t ds. 



Since 7 is a critical curve for X2 J t ds + X3 J ds among curves of length L, then we see that S £ = ii the 
ratio a : 2j3m is the same as A3 : A2. 

If 7 has inflection points, £ and the left hand side of (|2^) are still continuous functions on the space of 
rods, but J T ds can jump by ±7r under smooth variations of the curve (see [ pT| for a simple example). Thus, 
it only makes sense to prescribe A-tJj modulo tt. The gradient d J t ds still makes sense, and in fact the above 
proof goes through with this simple modification. □ 

We will now outline the computation of the first and second variation of energy £. Under our boundary 
conditions, the energy reduces to 

£ = an^ + pm^ ds, 



where m is the constant determined by (p4|). 

Let W he a periodic vector field along a closed curve 7, such that W preserves the length of the curve 
and respects the boundary conditions. Since any length-preserving deformation can be reparametrized so as 
to preserve arclength locally along the curve, we can assume, for the sake of convenience, that W preserves 
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a unit speed parametrization. Using variation formulas available in, for example, §6 of we see that W 
must be orthogonal to T, and the M^-derivative of m is given by 

Then 

£w = J < W", 2aKN > +2/3m < W' , kB > ds 
= J <W\-J > ds = J < W, J' > ds, 

where 

J = gT + 2aK'N + (2aKr - 2l3mK)B. 

The tangential component g of J can be chosen so that J' is normal to the curve and itself defines a 
length-preserving deformation of the curve, i.e., J < J",T > ds — Q. This gives 

g — an^ — A, 

where A is the appropriate constant. Then, by the usual reasoning, J' = for critical curves. 
Since the arclength derivative and the W^-derivative commute, 

£ww ^ J <W, VwJ' > ds = J - <W\ VwJ > 
Further computations yield 



ds. 



Eww ^ j 2a< W\ W" > +2(3m{< W, B >< W" , N> - < W\ B >< W, N >) 



+ (A - 3aK^) < W, W > ds+^ 

1j 



< W, kB > ds 



Our strategy for investigating stable non-planar rods will be to find conditions under which planar elastic 
rods are unstable. Letting W = fiN + vB, and assuming the torsion of the center line is zero, 

-\ 2 



2B 

-I- (A - 3aK^){^i^ + ly^) ds + -- 

1j 



L 

VK ds 



Of course, W must be the derivative of a periodic vector field along a planar curve, so /i and v are not 
arbitrary. 

Proposition 6.1. An elastic rod with circular centerline is stable if |r7i| < \/3 , and unstable if |m| > 
r-2-Ka 

This result is not new; it was derived, in a different form, by Zajac jist . For the sake of completeness, we 
give our own proof. 



Proof. Since curvature is a constant given by k = Stt/L, then A = an^ . If we set ^ = oq -t- a„ cos(na;) -I- 

n=l 

oo 

6„ sin(nx) and v = cq + c„ cos(nx) -f (i„ sin(nx), where x — ks, then the restrictions on W' force cq = 

n=l 

and ai —bi — 0. Now, 

oo 

£ww = 27r ^ [aK{n^{al + bl) + [n^ - + dl)) + 2l3mn{bnCn - a„d„)] 

n=2 

The term in the summation is a sum of two quadratic forms, in (a„, o?„) and in (6„, c„), which both have the 
same trace and determinant. These will be indefinite if (Q;K)^(n^ — 1) < /3^m^, and definite otherwise. □ 
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7. Stability of the Figure Eight 

The other planar closed centerline besides the circle is the figure-eight, where n — Kocn{-^,p). The 

non-constancy of the curvature makes it much more difficult to determine under what conditions £ww can 
be negative. We will proceed this way: treating £ww ^ functional on i-pcriodic functions /i(s) and j^(s) 
satisfying 

ip^ + iy'^)ds = l, (25) 

we will attempt to minimize £ww- The relevant Euler-Lagrange equations are 

(A - C - OK^)^ - 2afi" - 2(3mv' = 

(26) 



(A - C - 2,aK^)v - 



2av" + 2Pmfi' + —— / vKds = 0. 
L Jo 



Here, C is a Lagrange multiplier. One easily computes that, if ii and v satisfy these equations and (^5|) as 
well, then £ww = C. 

The above equations ( p6| ) uncouple if we make the 

Assumption 7.1. Our figure-eight elastic rod is untwisted {m — 0). 

In terms of the independent variable t = kqs/{2p), ( p6|) becomes 



'ep^sn^t- l + 4p2 + _L_ ,y + ^lJlcni / ly cntdt. (28) 



-„QL 

Here, we have used L — SpK/ng and A = KQa(2p^ — l)/(2p^) for the figure eight. 
Recall now Lame's equation: 

= {n{n + 1)/ sn^ t-h)y, (29) 

where n is usually taken to be a positive integer. The eigenvalues — that is, values for h for which there 
exist 4iir-periodic solutions — come in two flavours: for the lowest 2n + 1 eigenvalues, there exists a unique 
periodic solution (up to multiple), and this is given by a Lame polynomial, a polynomial of degree n in 
elliptic functions cn, dn and sn. For the remaining eigenvalues, which will not be of interest to us, periodic 
solutions coexist p^ ; that is, there always exist two independent periodic solutions. (How to write down the 
latter solutions in terms of theta functions is described in the last few pages of p4).) 

Proposition 7.2. Suppose W = fiN + vB for a periodic vector field W , and satisfies (27). If C < 0, 
then /i = 0. If C = then /i is a constant multiple of k, and £ww depends only on v. 

Proof. (27) is a Lame equation for n ~ 1. According to Ince the solutions of ( p9| ) for the three lowest 
eigenvalues are 

h=p' 
h = 1 




If C < 0, then h < 1 and ^ = dnt up to multiple, but this is not possible. For, if W — aT + bN + cB then 

a' = bn, b' — II — OK, c' = V. (30) 

Using K = cnt, and applying variation of parameters, we find that a and b will be periodic if and only if 
l"^^ fisntdntdt = and J^^ /i(l — 2p^ sn^ t) dt = 0. For /i = dni, the former is true but not the latter. 

If C=0, II = qn for some constant q, and W — qT will be a variation vector field which has the same action 
on £ as W, with (W - qT)' = vB. □ 

The analysis of the homogeneous equation (28), for the binormal component v of W' , involves the Lame 
operator for n~2: 

C := {d/dtf - {Qp^finH- h). 
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2p^C 

— 2 — makes C the left-hand side of (28). The solutions of Cy ~ for the 



y = sn i dn i 
y = cntdnt 


h=l + Ap^ 
h=\+p^ 


y = 1 - {1 + p^ - ^1 - p^ + p^) sn^ t 





Our objective will now be to obtain the value of [3 /a for which (28) has a nontrivial periodic solution for 
C = 0, and to show that this value is precisely the boundary between stability and instability. 

Lemma 7.3. For values of h lying in an open interval that extends to — oo and contains the three lowest 
eigenvalues (but no other eigenvalues) of £, there is a unique value of H for which 

l-iK 

Lv ~ H cnt I i^cntdt, (31) 
Jo 

has a nonzero periodic solution, and this value _ff is a continuous, strictly increasing function of h. 

Proof. First, assume that h is not an eigenvalue of £. Then the equation 

Cv = cnt. (32) 

has a unique periodic solution, constructed by variation of parameters. Any solution of ( ^l|) will be a multiple 
of this, and so H is uniquely determined, in terms of the solution of (|3^), by 

\ -1 



r4K 



H 



itdt 



Furthermore, if w — dv/dh for the solution of (p2|), then Cw = —v. Using the fact that C is self-adjoint. 



d 

dh 



4K 



4K 



iK 



ly cut dt = w Lv dt ~ ~ v Cw dt ~ — v dt 
/o Jo Jo Jo 

If h is either the first or third eigenvalue, ( p2[ ) has a solution, unique up to addition of an eigenfunction 
y, and H is determined as before. Moreover, when h is in the vicinity of either of these eigenvalues, we can 
construct a family of solutions to (32) depending continuously on h by requiring v to be an odd function of 
t — K . (This condition is already met by the non-eigenvalue solutions.) Thus, H depends continuously on h. 

li h = 1 no periodic solution to ( ^2| ) exists. Instead, in the vicinity of this eigenvalue, we let {v, H) 

be the unique solution to Cv ~ H cut such that v is periodic and J^^ ucntdt = 1. Then H depends 
continuously on h, with H = when h = 1 + p^ . □ 

Theorem 6. The untwisted figure eight elastic rod is stable if /J/a > 2 and unstable if /3/a < 2. 

Proof. We know, so far, that variations for which £ww < exist only if there is a solution of (|l]) for 

p'^fj . . ... . p'^P . 

Since iJ is a monotone increasing function of /i, this happens only when is 

Ka Ka 

1 + 4p^. We will now calculate that value. 



h <\ + Ap^ and H 



less than the critical value of H corresponding to h 

Let = sntdnt. Using reduction of order, we obtain a second, linearly independent solution of the Lame 
equation Cy = for h ~ 1 Ap^ , given by 



yi = cn^ tdnt 



P 



{p'Y 



cn t sn t — 



t + 



2p2 



1 



{p' 



n2 



dnt snt. 



where E{t) is an incomplete elliptic integral of the second kind Q]. 

From now on, let the prime denote d/dt. Rewrite the nonhomogenous equation 

dz 1' ' ' 

sn^t-h ; 



|) as the system 





r). 






\cntj ' 









Let Y = 



yi y2 



Vi 



y'2 



dt 



then variation of parameters gives z = Y [ £ + 



(33) 



where 



m 



Y-\x) 





cna; 



dx. 
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To ensure that z is periodic, we must take 

_ /p ^ yx cn t dt _ 

(Computations for the numerator and denominator both use the fact that, for the eUiptic modulus p associ- 
ated to the figure eight elastic rod, the complete elliptic integrals satisfy 2E — K.) Since the value of H is 
unaffected by adding multiples of 2/2 to ly, we can let C2 = 0. Then 

ycntdt= I {{ii{t) + ci)yi{t)+l2{t)y2{t))cntdt 

4K i-x ^4K 

ciitciix{y2(t)yi{x) - yi{t)y2{x))dxdt + ci / yi(t)cntdt 



Jo 







4K 

2 



?;i(t)(sn t + ci) cn tdt, 







using the fact that J y2 cnt dt — i sn^ t. Finally, we compute that v cn tdt = K/ {2p^), H = 2p^ / K, and 
/3/a = 2. □ 

Theorem^ implies that a conjecture made in Q and proved in Q, to the effect that there are no nonplanar 
stable closed elastic curves, would be false if extended to elastic rods. For, the twisted circle, with At/j — 2tt 
and m ~ 211 j L, is connected to the untwisted figure eight rod by a continuous 1-parameter family of rods. 
(A?/' — 2tt throughout the family up until we get to the figure eight, when Aip jumps to zero as the inflection 



points appear.) Proposition 6.1 and Theorem together imply that when \/3 < — < 2, both the circle and 

— . . . 

the figure eight rods are unstable for their respective boundary conditions. Since these two are the only 

closed planar elastic rod centerlines, it follows from the existence of a minimum that when (3 /a is in the 
above range there exist nonplanar closed elastic rods that are local minima for the elastic energy subject to 
the boundary condition = 2tt. 

References 

1. P. Byrd, M. Friedman, Handbook of Elliptic Integrals for Engineers and Physicists, Springer, 1953. 

2. A. Calini, T. Ivey, Bdcklund transformations and knots of constant torsion, to appear in J. Knot Theory and its Ramifi- 
cations (1998). 

3. W. Ince, The Periodic Lame Functions, Proc. Royal Soc. Edin. 60 (1940) 47-63. 

4. S. Kida, A vortex filament moving vuithout change of form, J. Fluid Mech 112 (1981), 397-409. 

5. J. Langer, D. Singer, Knotted Elastic Curves in R^^ J. London Math. Soc. (2), 30 (1984), 512-520. 

6. — , Lagrangian Aspects of the Kirchhoff Elastic Rod, SIAM Review 38 (1996), 605-618. 

7. — ,Curve Straightening and a Minimax Argument for Closed Elastic Curves, Topology 24 (1985), 75-88. 

8. Y. Li, J. Maddocks, On the Computation of Equilibria of Elastic Rods; Part I: Integrals, Symmetry and a Hamiltonian 
Formulation, preprint. University of Maryland (1996). 

9. R. Manning, K. Rogers, J. Maddocks, Isoperimetric Conjugate Points uiith Application to the Stability of DNA Minicircles, 
preprint. University of Maryland (1997). 

10. D. Mumford, Elastica and Computer Vision, pp. 491—506 in "Algebraic geometry and its applications", ed. C. L. Bajaj, 
Springer, 1994. 

11. R,. Ricca, H. Moffatt, The Helicity of a Knotted Vortex Filament, pp. 225—236 in "Topological Aspects of the Dynamics of 
Fluids and Plasmas", ed. Moffatt, Kluwer, 1992. 

12. K. A. Rogers, Stability Exchange in Parameter- Dependent Constrained Variational Principles with Applications to Elastic 
Rod Models of DNA Minicircles (thesis). University of Maryland, 1997. 

13. W. Magnus, S. Winkler, Hill's equation, Interscience, 1966. 

14. E.T. Whittaker, G.N. Watson, A Course of Modern Analysis, Cambridge, 1927. 

15. E. E. Zajac, Stability of Two Planar Loop Elasticas, Transactionsn of ASME Series E, Journal of Applied Mechanics 29 
(1962), 136-142. 

(Ivey) Dept. of Mathematical Sciences, Ball State University, Muncie IN 47306 
E-mail address: tivey@math.bsu.edu 



(Singer) Dept. of Mathematics, Case Western Reserve University, Cleveland, OH 44106-7058 
E-mail address: das5Spo.cwru.edu 



